CMB Anisotropy in Compact Hyperbolic universes I: Computing Correlation 

Functions 



On 



J.Richard Bond, Dmitry Pogosyan and Tarun Souradeep 
Canadian Institute for Theoretical Astrophysics, 
University of Toronto, ON M5S 3H8, Canada 
(February 1, 2008) 

CMB anisotropy measurements have brought the issue of global topology of the universe from 
the realm of theoretical possibility to within the grasp of observations. The global topology of the 
universe modifies the correlation properties of cosmic fields. In particular, strong correlations are 
predicted in CMB anisotropy patterns on the largest observable scales if the size of the Universe is 
comparable to the distance to the CMB last scattering surface. We describe in detail our completely 
general scheme using a regularized method of images for calculating such correlation functions in 
models with nontrivial topology, and apply it to the computationally challenging compact hyperbolic 
spaces . Our procedure directly sums over images within a specified radius, ideally many times 
the diameter of the space, effectively treats more distant images in a continuous approximation, 
and uses Cesaro resummation to further sharpen the results. At all levels of approximation the 
O ■ symmetries of the space are preserved in the correlation function. This new technique eliminates 

the need for the difficult task of spatial eigenmode decomposition on these spaces. Although the 
U — i , eigenspectrum can be obtained by this method if desired, at a given level of approximation the 

^ ■ correlation functions are more accurately determined. We use the 3-torus example to demonstrate 

that the method works very well. We apply it to power spectrum as well as correlation function 
evaluations in a number of compact hyperbolic (CH) spaces. Application to the computation of 
' CMB anisotropy correlations on CH spaces, and the observational constraints following from them, 

are given in a companion paper. 

(N ' 

^ i. introduction 

On 

0^ ■ The remarkable degree of isotropy of the cosmic microwave background (CMB) points to homogeneous and isotropic 
Freidmann-Robertson- Walker (FRW) models for the universe. The underlying Einstein's equations of gravitation are 
purely local, completely unaffected by the global topological structure of space-time. In fact, in the absence of spatially 
inhomogeneous perturbations, a FRW model predicts an isotropic CMB regardless of the global topology. 

Flat or 'open' FRW models adequately describe the observed average local properties of our Universe. Much of 
recent astrophysical data suggest the cosmological density parameter in nonrelativistic matter, f2 m , is subcritical |3|. 
The total density parameter f2o includes relativistic particles and vacuum, scalar field or cosmological constant, 
■ contributions, as well as fl m . If a cosmological constant (or some other exotic smooth matter component) does not 
, ' compensate for the deficit from unity, this would imply a hyperbolic spatial geometry (uniform negative curvature) , 
commonly referred to as an 'open' universe in the cosmological literature. The simply connected (topologically trivial) 
hyperbolic 3-space Ti. 3 and the flat Euclidean 3-space £ 3 are non-compact and have infinite volume. There are 
numerous theoretical motivations, however, to favor a spatially compact universe To reconcile a compact 

universe with a flat or hyperbolic geometry, consideration of multiply connected (topologically non-trivial) spaces is 
required. (Inhomogeneous simply connected models are another way out; e.g., a hyperbolic bubble could be embedded 
in a highly inhomogeneous space and if the bubble is much larger than the Hubble radius we would not know.) Compact 
hyperbolic spaces have been recently used to construct cosmological models within the framework of string theory |j) . 
A compilation of recent articles describing various aspects of current research in compact cosmologies can be found 
in@. 

The realization that the universe with the same local geometry has many different choices of global topology is 
as old as modern cosmology - De Sitter was quick to point out that Einstein's closed static universe model with 
spherical geometry S 3 could equally well correspond to a multiply connected Elliptical universe model where the 
antipodal points of S 3 are topologically identified. The possibility of a multiply connected universe has lingered 
on the fringes of cosmology largely as a theoretical curiosity (e.g., in Q), but was sometimes invoked to explain 
puzzling cosmological observations (e.g., as a possible explanation for the isotropy of the CMB radiation [Q, and for 
the (controversial) observations of periodic or discordant quasar redshifts Jll|]). There is a long history of attempts 
to search for signatures of non-trivial global topology by identifying ghost images of local galaxies, clusters or quasars 
at higher redshifts [jl,^ 12 1. The search for signatures of global topology in the distribution of luminous matter can 



probe the topology of the universe only on scales substantially smaller than the apparent radius of the observable 
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universe. Another avenue in the search for global topology of the universe is through the effect on the power spectrum 
of cosmic density perturbation fields, reflected in observables such as the distribution of matter in the universe and 
the CMB anisotropy. 

The observed large scale structure in the universe implies spatially inhomogeneous primordial perturbations existed 
which gave rise to the observed anisotropy of the CMB. The global topology of the universe also modifies the local 
observable properties of the CMB anisotropy on length scales up to a few times the horizon size. In compact universe 
models, the finite spatial size usually precludes the existence of primordial fluctuations with wavelengths above a 
characteristic scale related to the size. As a result, the power in CMB anisotropy is suppressed on large angular 
scales. Another consequence is the breaking of statistical isotropy in characteristic patterns determined by the photon 
geodesic structure of the manifold. One can search for such patterns statistically in CMB anisotropy maps. Full-sky 
CMB data, such as from the COBE-DMR experiment, can constrain the size of the universe and its topology to the 
extent that such correlation patterns are absent in the data |2|,[l6| . 

For Gaussian perturbations, the angular correlation function, C(q,q'), of the CMB temperature fluctuations in 
two directions q and q' in the sky completely encodes the CMB anisotropy predictions of a model. For adiabatic 
perturbations, the dominant contribution to the anisotropy in the CMB temperature measured with a wide-angle beam 

1/2 

(#fwhm ^ 2° fi ) comes from the cosmological metric perturbations through the Sachs- Wolfe effect. The angular 
correlation function C(q,q r ) then depends on the spatial two point correlation function, = ($(x, t ls )<I>(x', t ls )), of 
the gravitational potential, on the three- hypersurface of last scattering along the lines-of-sight, q and q' . Thus we 
need to learn how to compute spatial correlation functions on compact spaces. 

When the eigenfunctions of the Laplacian on the space are known, the correlation function can be readily obtained 
via a mode sum. For this to be tractable the known eigenfunctions would, preferably, be expressed in a reasonably 
simple closed form. However, obtaining closed form expressions for eigenfunctions may not be possible beyond the 
simplest topologies. Some examples where explicit eigenfunctions have been used include flat models |l7]-|2(J and 
a non-compact hyperbolic space with a horn topology In this paper we provide a detailed description of our 

regularized method of images for computing correlation functions, which does not require any prior knowledge of the 
eigenfunctions on the compact spaces j^|l| . The method allows us to accurately compute the correlation function on 
compact spaces where eigenfunctions are not available. Most notorious in this respect are compact spaces with uniform 
negative curvature, the Compact Hyperbolic (CH) spaces, for which even numerical estimation of the eigenfunctions is 
believed to be a challenging task (for recent progress see |13-15|). A novel feature of our scheme is the regularization 
procedure that we devise in order to successfully implement the method of images for the power spectra of cosmological 
perturbations expected from early universe physics. As an additional bonus, this regularization scheme enhances the 
convergence properties of the method, which proves to be very useful for tackling CH spaces. 

In §fll, we introduce compact spaces and briefly review the aspects that are relevant for our work. The method of 



images is derived in § [II . As a simple illustrative example, we apply the method of images to the case of a simple flat 
torus in § [V . The implementation of the method in CH spaces is presented in §0. T he gr oss properties of the power 
spectrum in CH spaces that can be inferred from our computation is discussed in §VC. We derive the properties 



of the correlation function in CH space in § VI. In § VII , we discuss our results and point to our companion paper 



describing our detailed CMB anisotropy calculations fllf 



II. PRIMER ON COMPACT SPACES 



A. Mathematical Preliminaries 



A compact cosmological model M. is a Quotient Space constructed by identifying points on the standard FRW 
spaces under the action of a suitable [] discrete subgroup of motions, T, of the full isometry group G, of the FRW 
space (see p|p^]). The isometry group G is the group of motions which preserves the distances between points. The 
infinite FRW spatial hypersurface is the universal cover, M u , tiled by copies of the compact space, M.. Cosmological 
postulates of local homogeneity and isotropy restrict M. u to be a space of constant curvature - Hyperbolic TC 3 , with 
negative curvature, Spherical S 3 , with positive curvature, or flat Euclidean £ 3 , with zero curvature. The compact 



1 A Quotient space M = M u /V consists of the equivalence classes of points of M u equivalent under F. M is a manifold 
when r is a subgroup of the isometry group, G of A4 U which acts freely (fixed point free: 7X = x 7 = J) and properly 
discontinously (every point x £ M has a neighborhood t/ x such that 7(t/ x ) (7 X = 0V7 £ F, 7 / /). The manifold Ai is 
compact if the corresponding Dirichlet domain is compact. 
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space for a given location of the observer is represented as the Dirichlet domain with the observer at its basepoint. 
Every point x of the compact space has an image Xj = 7^x in each copy of the Dirichlet domain on the universal 
cover, where ji £ Y. The tiling of the universal cover with Dirichlet domains is a Voronoi tessellation (a familiar 
concept in cosmology that has been used in modeling the large scale structure in the universe) , with the seeds being 
the basepoint and its images. The Dirichlet domain represents the compact space as a convex polyhedron with an 
even number of faces, with congruent pairs of faces identified (glued) under Y. For more details, see e.g., |i|, |2^ , |23| ,|6f . 

More explicitly, the Dirichlet domain around the basepoint, xo, is the set of all points on the universal cover which 
are closer (or equidistant) to Xo than to any of the images 7x0 (7 £ Y) of the basepoint. The hyperplane that bisects 
the segment joining x to 7x0 divides M u into two parts; let H 7 denote that half that contains the basepoint. By 
definition, the Dirichlet domain around Xo is 2? xo = C\^H 1 . Thus, the Dirichlet domain is bounded by hyperplanes 
bisecting the (geodesic) segments joining the basepoint x to a set of adjacent images. (The corresponding set of 
motions is called the set of adjacency transformations.) The faces of the polyhedron are identified pairwise; the one 
formed by the hyperplane bisecting the segment joining xo to 7x0 is identified with the corresponding one bisecting 
the segment joining xo to 7~ 1 xo- The adjacency transformations generate the group Y and hence are also known as 
the face generators. 

In the context of cosmology, the Dirichlet domain constructed around the observer represents the universe as 'seen' 
by the observer. (In our papers, we will often loosely use the same notation M. to refer to the compact space as well 
as one of its Dirichlet domain representations whose basepoint is clear from the context.) It proves useful then to 
define the outradius R > and the inradius R < of the Dirichlet domain j^]: the radius of the smallest sphere around 
the observer which encloses the Dirichlet domain and that of the largest sphere around the observer which can be 
enclosed within the Dirichlet domain, respectively . The ratio i?> / i?< is a good indicator of the shape of the Dirichlet 
domain. 

In cosmology, distances are inferred from the cosmological redshift which is related to the light travel time. Con- 
sequently, in compact spaces, and, more generally, in multiply connected spaces, light from the same source seems to 
arrive from points at different locations - the source and its images on the universal cover. Any source inferred to be 
at a distance greater than R > from the observer is bound to be the image of a point which is physically closer. On 
the other hand, a source closer than i?< is definitely at its true physical distance. Note that i?> and i?< are specific 
to the location of the observer within the compact space since the Dirichlet domain around different observers can, 
in general, vary. An observer (and Dirichlet domain) independent linear measure of the size of the compact space is 
given by the diameter of the space, o\m = sup^g^ d(x, y), i.e., the maximum separation between two points in the 
compact space. 

The isometry group G defining the global symmetries of M is the centralizer of Y in the isometry group G u of 
its universal cover, M u ; i.e., G — {g £ G u \g^ — 73V7 € Y}. In general, M respects less symmetries than A4 V ; 
consequently, G is of lower dimension than G u . M. is (globally) homogeneous if and only if G is transitive on Ad: 
i.e., for any two points x, y £ M., there exists a g £ G such that y = gx (an equivalent statement is that a compact 
space is globally homogeneous if and only if every element of 7 £ Y is a Clifford translation, i.e., the displacement 
function, d 7 (x) = d(x, 7X), is independent of x for all points x £ M. u ). Space is isotropic (around a point, or observer, 
x ) if G contains a subgroup of rotations around x Q . The only example of a multiply connected compact universe 
which retains all the symmetries of its universal cover is the Elliptic space (S 3 /Z2) with spherical geometry; the simple 
flat torus is anisotropic, and all others (including the entire class of CH manifolds) break global homogeneity as well. 
As the result, the two-point correlation functions £(x, y) in CH spaces depend separately on both points x and y, 
and not only on the distance c?(x,y,) as in the familiar FRW spaces. 

B. Compact Euclidean spaces 

The compact spaces with Euclidean geometry (zero curvature) have been completely classified. In three dimensions, 
there are known to be six possible topologies that lead to orientable spaces p], p2| , p3"| : the simple flat torus where £ 3 is 
identified under a discrete group of translations; three other flat tori where the identification is under a screw motion, 
i.e., translation accompanied by a rotation about the direction of the translation, namely, rotation by tt or tt/2 in 
one component of the translations, and rotation by tt in all three directions; and finally, two topologies made by 
identifying the planar hexagonal lattice under a screw motion perpendicular to the plane with rotation of 7r/3 and 
27r/3, respectively. 
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C. Compact Hyperbolic spaces 



For cosmological CH models, M. u = Ti 3 , the three dimensional hyperbolic (uniform negative curvature) manifold 
with line element 

ds 2 = d X 2 + sinh 2 x (d9 2 + sin 2 6d(f> 2 ) , (1) 

where \ — ( T o — T)/d c is the affine distance, r is the conformal time and d c — cHq 1 A/l — f^o is the curvature radius 
set by the present cosmological density parameter fio and the Hubble constant Hq.h For a universe with hyperbolic 
geometry, < Oq < 1; thus the size of d c ranges from cHq 1 as f^o — > to infinity as flo — > 1. 

7i 3 can be viewed as a hyperbolic section embedded in four dimensional flat Lorentzian (Minkowski) space, by 
representing each point on Ti 3 as a unit four vector, X = (£(x), x), normalized by d c in Minkowski space (£ 2 — |x| 2 = 
d 2 ). The distance between two points on Ti 3 is given by the dot product of the corresponding four vectors. The 
isometry group of TL 3 is then the group of rotations in the four space, the proper Lorentz group SO(3, 1). A CH 
manifold is completely described by a discrete subgroup, T, of the proper Lorentz group 30(3, 1). 

There are two remarkable features of tessellating Ti 3 under a discrete group of motions that are absent in the 
flat geometry. Firstly, whereas in flat geometry all finite volume Quotient Spaces obtained are necessarily compact, 
one can tessellate Ti. 3 with tiles of finite volume which are non-compact, giving rise to a class of non-compact finite- 
volume hyperbolic universe models. Typically, these spaces have cusp- like extensions to infinity. Secondly, a given 
CH topological structure fits only for a specific volume of the space. This is in contrast to flat compact spaces, in 
which the same topological structure can be imposed on any scales. In particular, all simple flat tori with different 
identification lengths are homeomorphic to each other, i.e., one can be obtained from the other by a continuous 
mapping/deformation, but are not isometric to each other, i.e., the distance between the mapped points are not the 
same. Two hyperbolic spaces of finite volume which are homeomorphic are necessarily isometric (in three dimensions 
and above). This result, known in mathematics literature as the strong rigidity theorem, can be attributed to the 
existence of the intrinsic length scale on TL 3 - the curvature radius, d c . It not difficult to convince oneself that the 
Dirichlet domain has to have a fixed size relative to d c in order to tile Ti 3 - if one were to scale the polygon size 
(relative to d c ) then the scaled tiles will no longer fit together at the vertices. 

Thus, a CH manifold, M, is characterized by a dimensionless number, Vm = VM/d 3 , where Vm is the volume of 
the space and d c is the curvature radius f24j] . There are a countably infinite number of CH manifolds, and no upper 
bound on V M . The theoretical lower bound stands at V M > 5/(2V3)Arcsinh 2 (\/3/5) « 0.167 [{25). Tne smallest CH 
manifold discovered so far has Vm = 0.94 and it has been conjectured that this is, in fact, the smallest possible p6"| , f27f . 
The physical volume of the CH universe with a given topology, i.e., a fixed value of Vjn/d 3 ,, is then set by d c and is 
thus related to cosmological parameters. The Geometry Center at the University of Minnesota has a large census of 
CH manifolds in its public domain. This census was created using the SnapPea computer software which is also freely 
available at the website p8| . The Minnesota census lists several thousands of these manifolds with Vm up to ^ 7 and 
the SnapPea software can be used to obtain various characteristic properties of these CH manifolds such as the ones 
listed in Table | and also the generators of the discrete group motion T that we need for our computational method. 

The sequences of CH manifolds are closed well-ordered sets of order type w", i.e., are arranged in a countably 
infinite number of countably infinite sequences (in increasing volume), with each sequence having a non-compact 
finite volume hyperbolic space as its limit. These limit spaces of finite volume have cusps extending to infinity. The 
cusps are homeomorphic to T 2 x R + , the product of a 2-torus and the positive real line. The sequence of CH spaces 
arises through Dehn filling, a procedure which truncates the cusps of the limit space along a closed curve on T 2 (with 
allowed pair of winding numbers) and glues in a solid torus (by identifying the closed curve along a meridian of the 
solid torus). The nomenclature of CH spaces in the Minnesota census, which we follow in this work, identifies a 
sequence by the limit space and specifies a CH space belonging to it by two (coprime) integer indices in parenthesis 
correspond to the winding numbers, to and n, which characterize the closed curve on the torus; e.g., m004(-5,l) refers 
to the CH manifold obtained by Dehn-filling the single cusp of the non-compact space m004 by a torus glued along a 
closed curve with winding numbers (—5, 1). The limit non-compact space is achieved as a limit of large values of the 
winding numbers in the corresponding sequence of CH manifolds. From the observational perspective, the CH spaces 
with high winding numbers should be practically indistinguishable from the limit non-compact space. The space with 
winding numbers (to, n) is equivalent to the space (—to, —n) in the same sequence. Not all integer values of m and n 
lead to a CH space, there are bounded forbidden gap regions in the m-n lattice. One should note that the volume 
does not uniquely characterize a space; there can be a finite number of distinct CH spaces with the same volume 



2 Unless explicitly indicated, all distances and times in non-flat geometry are in units of the curvature scale d c . 
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(besides the mirror image pairs for chiral CH spaces such as m004(-5, 1) & m004(5, 1)). Also, the nomenclature 
described above is not unique: the same CH space may belong to two different sequences, e.g., the space m003(-2,3) 
is equivalent to m004(5,l). (See |^7j for details.) 

All CH hyperbolic spaces are necessarily globally inhomogeneous since the only element of 50(3, 1) which is a 
Clifford translation is the identity p2[ . A short proof follows. Assume that g <E 50(3, f) is a Clifford translation. 
Then d(xi,<7Xi) = d(x2,(?X2) for any two points xi and X2 on 7i 3 . Since 5*0(3,1) is the isometry group of the 
homogeneous space H 3 , there exists a Lorentz transformation A 6 50(3,1) such that X2 = Axi. Consequently, 
d(xi,gxi) = d(Axi, gAxi), implying that g = A~\gA VA e 50(3,1). Since the Lorentz group is non-Abelian 
(i.e., actions of the group elements do not commute), g must be the identity. 




FIG. 1. The Dirichlet domains of two of the compact hyperbolic (CH) spaces that we have studied are shown. On the left 
is a small CH space, m004(-5,l), and on the right is a large CH space, v3543(2,3). 



TABLE I. The characteristics of some of the compact hyperbolic (CH) manifolds that we have studied. The nomenclature 
of the spaces conforms to the Minnesota census. The volume and diameter (expressed in the units of curvature radius) are 
topological invariants which relate to the density of states in those spaces. The specifications have been obtained from the 
Minnesota census using SnapPea. The exception is the diameter which is our estimate using a simple random sampling 
algorithm. However, in the context of cosmology, the observer specific characterization of the Dirichlet domain (around the 
observer) proves useful. The m003 models correspond to spaces obtained by different Dehn fillings on a (non-compact) cusped 
manifold, m003. The CH space m003(-3,l) has the smallest volume in the m003 series and is also currently the smallest CH 
space known. The m004(-5,l) has the smallest volume in the m004 series. The v3543(2,3) space is a relatively large one. 



Properties 


m003(-3,l) 


m003(-4,3) 


m003(-4,l) 


m003(2,3) 


m003(-5,4) 


m004(-5,l) 


v3543(2,3) 


Topological invariants: 
















Volume: Vm 


0.94 


1.26 


1.42 


1.54 


1.59 


0.98 


6.45 


diameter: d,M 


0.84 


1.01 


1.10 


1.16 


1.21 


0.86 


1.90 


First Homology group 


z 5 ®z 5 


Z 5 ®Z 5 


Z35 


Z35 


Z30 


z 5 


Z93 


Dirichlet domain specific: 
















outradius: R > /d c 


0.75 


0.84 


1.07 


0.83 


0.94 


0.75 


1.33 


inradius: R < /d c 


0.52 


0.55 


0.54 


0.59 


0.57 


0.53 


0.89 


No. of faces 


18 


22 


24 


26 


28 


16 


38 


No. of vertices 


26 


40 


44 


48 


52 


16 


72 
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III. METHOD OF IMAGES 



A. The correlation function 



The correlation function £$(x, x') of a scalar field, can be expressed formally in terms of the orthonormal set of 
eigenf unctions ^ of the Laplace operator, V 2 , on the hypersurface (with positive eigenvalues kf > 0), as p9| 

e*(x,x') ^ftWp^x]^), where (V 2 +fc?)% =0. (2) 

» 3=1 

The spectrum of the Laplacian on a compact space (thus with closed boundary conditions) is a discrete ordered set of 
eigenvalues {kf} (fc 2 = and kf < kf +1 ) with multiplicities m,. The function P$(fcj) describes the rms amplitude of 
the eigenmode expansion of the held <£, determined in the context of cosmology by the physical mechanism responsible 
for the generation of 

Except in simple cases, neither the spectrum of the Laplacian {kf} nor the eigenfunctions \I/£-(x) are known for 
compact manifolds, so eq.(0) cannot be used to calculate the correlation function £J(x, x') on a compact manifold 
M. directly. In contrast, the universal cover M. u for the compact manifold is usually simple enough (e.g.,Tt 3 , S 3 or 
£ 3 ) that the eigenfunctions \l/"(fc,x) are known and the correlation functions £$(x, x') are easily computable. We 
consider H 3 and £ 3 geometries where the spectrum of eigenvalues on the universal cover M. u is continuous, reflected 
in the notation \&"(fc,x) replacing the discrete index % with a functional dependence on k. 

The regularized method of images we developed in allows computation of the correlation function £J(x, x') on 
a compact (more generally, a multiply connected) manifold M. — A4 U /T from the correlation function, £|(x, x'), 
on M u . We now give the explicit derivation of the relation between £J(x,x') and £|(x, x'), calculated with the 
same form of the power spectrum P$(fc). The method solely relies on the knowledge of the action of elements of the 
discrete group, T, and requires no information regarding the eigenvalues and eigenmodes of the Laplacian on M. The 
expansion (||) is equivalent to the system of integral equations on the correlation functions, £| and 

/ dx'&(x,x') *f i (x')=P$(fc i )^ i (x), (3a) 

/ dx' e$(x,x') * 3 "(fc,x') = P,(fc)tf?(fc,x) . (3b) 

Eq.(|t]) is a consequence of the homogeneity of .M" through a theorem j2{| which states that the eigenfunctions of 
the Laplacian are also the eigenfunctions of the integral operator corresponding to any two-point function which is 
point-pair invariant, i.e., depends only the distance between the points. The orthonormality of the eigenfunctions 
leads to the expansion (|3a|) for 

Every eigenfunction vf^- (x) on M. is also an eigenfunction of the Laplacian on the universal cover A4 U with eigenvalue 
kf, hence is a linear combination of degenerate eigenfunctions ^"(fc^x) on M. u with eigenvalue kf, i.e., \&f,(x) = 
Si' a i,ji'^'(^») x )' Thus a subset of equations ( pb| ) can be written as 

/ dx' £«(x,x') ^.(x') - P 4 (k)*£(x) ( 4 ) 

= / dx'£|(x,x') ^.(x'). 
Using automorphism of with respect to T, vf^-^x) = \t?-(x)V7 6 L, and the fact that M. tcsscllates M. u , 



f dx' £J(x,x') ^.(x') = / dx' a(x,x')^(x') = W dx' &(x, 7 x') ^(x') (5a) 

^,(x'), (5b) 



dx' 
M 



7er 



where ^2 denotes a possible need for regularization at the last step when the order of integration and summation is 
reversed. Since cq.(5b) is satisfied for all comparison with eq. ( |3a| ) implies that 



G 



a(x,x') 



7er 



,7x')- 



(6) 



This is the main equation of our method of images, which expresses the correlation function on a compact space 
(and more generally, any non-simply connected space) as a sum over the correlation function on its universal cover 
calculated between x and the images jx' (7 S T) of x'. 

Regularization in eq.(5b) is required when the correlation function on the universal cover does not have a compact 
support, J MU £|(x, x')c?x' = 00. Let us consider a new two-point function on M u defined by 



a(x,xo = a(x,x')-^ 



dx" a(x,x") 



(7) 



7A-1 



for 7 such that x' lies in 7 At. Replacing £|(x, x ') by its regularized version £$(x, x') in eq. (|5a|) gives the same 
integral equation on £|,(x, x') as the original one for every eigenfunction vE^- with the exception of that for the 
zero mode k 2 — 0. For the zero mode, = const, for which the right-hand side of eq. (0) is now zero 
Im" C$( x , x')c?x' = 0. Thus, the regularized sum over images is 



Also, 



a (x, x') = x: a (x, tx') = J2 ^ ( x > t x ') - 77— / dx" a (x, x") 

7 er 7G r Vm Jm " 



(8) 



One can verify that £S( X , x') is biautomorphic with respect to T, i.e., £|,(x, x') = £|(7ix, 72 x ') V7i,72 G T. In 
addition, £| is smooth and symmetric, £|,(x, x') = £|,(x',x), since £| is smooth and point-pair invariant. With these 
conditions satisfied, the correlation function then qualifies to be the kernel of an integral operator on square integrable 
functions on M.. 

The sum in eq. (|^) can be obtained analytically in a closed form only in a few cases, e.g., the simple flat torus (see 
section IV). For a numerical implementation of the sum over images, it is useful to present the formal expression (|^) 
as the limit of a sequence of partial sums 



£S(x,x') = lim 

N— J-O 



N 

X 

i=0 



#(x, 7i x')--i- / dx"8(x, 7 ,x") 
Km Jm 



(9) 



with the discrete motions 7, sorted in increasing separation, <i(x, 7^x') < d(x, 7j +1 x'). In practice, the summation over 
images is carried out over a sufficiently large but finite number of images, from which the limit N — > 00 is estimated. 
Accuracy is enhanced if the ji used correspond to a tessellation of A4 U with the basepoint of the Dirichlet domain 
shifted to x. 

A more effective and simpler limiting procedure for implementing the regularized method of images is to explicitly 
sum images up to a radius r* and regularize by subtracting the integral £|(r) over a spherical ball of radius r». 
This further eliminates the need to know the precise shape of Dirichlet domains (and integrate o ver potentially 
complicated shapes). We use this limiting procedure when dealing with CH spaces (see eq. ( p4| ) in § VB). The example 
of the simple flat torus discussed in §IV shows that the radial limiting procedure does better in eliminating unwanted 
low k power than eq. (^). On the other hand, the radial limiting procedure cannot be applied unambiguously if the 
number of images used is small whereas eq.(||) holds at all values of N. 

As we see, the need for regularization is strictly dictated by the form of £ u , which encodes how modes are excited, 
and is not specific to CH spaces. Some authors |5C| have incorrectly attributed the need for the regularization that 
we invoke in [[110 to the exponential proliferation of periodic orbits or the chaotic nature of classical trajectories in 
CH spaces. Regularization is required in flat compact spaces too (see §IV) if the £ u mode expansion contains the zero 
mode k 2 — 0. On the other hand, if has compact support there is no need for regularization even in CH space; 
this, in fact, holds under weaker conditions. 

The counter term in eq.(^J) significantly improves the convergence of the sum over images, even if in the limit 
N — > 00 the term is zero ( i.e., J Mu f$( x , x')<ix' = 0) and regularization is formally not required. Indeed, in this case 
the A^th partial term of eq.(tt) is equal to 



3 The volume integral of the Laplace equation on a compact manifold M (thus with closed boundary) gives kl J M dx'^^(x) 
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a(x,x') k =^a(x,7ix') + t- / dx" a(x,x") . (10) 

where M.m — A4 U — ufL {~fiA4} is the complement relative to M. u of the domains from which the image contribution 
has been explicitly summed. Thus, £$(x, x')| N corresponds to the approximation where the first N images up to 
the distance d < d(x, 7iv x ') are summed explicitly and the contribution of the rest of the images is estimated as 
the integral. The latter estimation is quite natural, since the sum over densely packed distant images is similar to a 
Monte-Carlo expression for the integral. 

Finally, we note that the regularization term is certainly not unique. In fact, instead of starting with a regularized 
correlation function on the universal cover as in eq. (Q), we could have started with a regularized scalar field, $ = 
$ — J M dV^/VM- This gives rises to different regularizing counterterms that might possibly be more effective but 
are also more complicated. In the same spirit of viewing the field as the primary starting point, the sum-over-images 
representation for £J(x, x ') can be viewed as a double sum £$(7 X , 7'x'). In general this is computationally 

more expensive. However, for the case of the simple flat torus, the double sum can be expressed as a single summation 
over appropriately weighted image contributions and has a significantly faster convergence. 



B. The Power Spectrum 



An illuminating way to present the power spectrum of fluctuations in compact spaces is to separate the density 
of states (the spectrum of eigenvalues weighted with multiplicity), solely determined by the geometry and topology 
of the space, from the occupation number (rms amplitudes of the eigenmodes) determined by the physics invoked to 
excite the available modes. This allows a discussion of the effect of non-trivial topology independent of the generation 
mechanism and the resultant statistical nature and spectrum of the initial perturbations. 

We define the power spectrum "P|,(fc,x) of the scalar field $ to be the function which describes the contribution 
to the variance of the field cr|,(x) from the modes in a logarithmic interval of eigenvalues d\nk. The variance is 
given by the correlation function at zero lag, thus our definition satisfies £$(x, x) = J Q d\n(k)T'^(k,x). The power 
spectrum depends both on the eigenvalue spectrum of the Laplacian, described through a collapsed two point function, 
n c (fc,x, x), as well as on the rms amplitudes P$(k) of the eigenmode expansion of the field on the universal cover: 

P£(fc,x)dln(fc) = P< s ,(k)n c (k,x,x)dln(k), (11a) 
n c (fc,x,x) = kY, 6{k- ki)S ki (x), (lib) 

i 

|2 



sm = J2\*U*)\ ■ ( llc ) 
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The function n c (fc,x, x) itself may be interpreted as the local density of states (per dhik) per unit volume at the 
point x. The notation n c (fc,x, x) explicitly retains the intimate connection to a two-point kernel. Indeed, Sfc ( (x) is 
related to the imaginary part of the Green function 

Ik 

S ki (x) = ^Im[G fc4 (x,x)]. (12) 

7T 

The dependence on the position is a manifestation of the global inhomogeneity of the space. In the case of a globally 
homogeneous space, Gfc 4 (x,x') = Gk f (x — x'), implying that 5fc 4 (x) and n(k,x) are position independent. We shall 
denote mean density of states by n c (k) 7 just omitting spatial dependence. Integrating Sfe^x) over the volume of the 
manifold M. gives the multiplicity of the ith eigenvalue J M 5fe i (x)dx = to.;, hence 

n c (k) = — !— / n c (fc,x, x) dx = — — k S(k — ki)rrii . (13) 
Vm Jm i 

Thus we factor the power spectrum in a compact space into a rms amplitude of modes and the density of states: 



n(*0 = T7-/ V%(k,x)dx^P^{k)n c {k). (14) 
>M 



In the universal covering space M u with a continuous spectrum of eigenstates, the density of states kernel is given 
by 
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n u (k, x, x') = — ^ £ *«(fc 7 x )*«* (fc, x ') . (15) 
i 

The homogeneity of A4 U implies that n"(fc,x, x) = n u (k) = fc 3 /(27r 2 ) is independent of x. Akin to the correlation 
function, the method of images can be used to establish the connection between the local density of states kernel on 
the compact manifold Ai c and on its universal covering space A4 U : 

n c (fc,x,x) = ^V(fc,x, 7 x). (16) 

7 

The method of images applied to the density of states, together with the connection to the trace of the Green 
function from eq.([l^), is what is essentially embodied in the celebrated Selberg trace formula [p9[ . We emphasize 
that, although the basic ideas are similar, the computation of the correlation function and the density of states n c (k) 
in a compact space are distinct problems. Getting the density of states is computationally more challenging since 
the aim is to recover a singular function (string of delta functions) by superposing smooth functions. For computing 
correlations of the CMB anisotropy in compact spaces, we only need to apply the method of images to the correlation 
function. Also, the pairwise correlation function at distinct points x ^ x', calculated by the method of images at any 
level of approximation, satisfies exactly the periodicity of the space, which the density of states does only if determined 
with absolute precision. 



IV. FLAT TORUS MODEL: A SIMPLE EXAMPLE 



The simple fiat torus model, T 3 , is the compactification of the 3-dimensional Euclidean space which identifies 
points under a discrete set of translations, x — * x + nL, where L is the size of the torus and n is a vector with integer 
components. The corresponding Dirichlet domain is a cube (more generally, a parallelepiped) with opposite faces 
identified (glued together). This is, in fact, the model one is studying when one simulates the universe in a finite box 
with periodic boundary conditions. 

Since the eigenfunctions of the Laplacian on T 3 are simply the discrete plane waves, the evaluation of £|(x, x ') on 
T 3 as a sum over modes functions (see eq.(||)) using Fast Fourier Transforms is a preferred technique, one we have 
used extensively to constrain the size of such models using the COBE-DMRdata. However, we revisit this simple case 
to illustrate the various steps and clarify subtleties involved in the calculation of £|(x, x') in a multiply connected 
universe using the method of images. The toroidal case also provides a good benchmark for evaluating the efficiency 
of the method of images. 

The correlation function in the periodic box implied by the T 3 topology is 



£f ( X . X ') = ^X! p <f( k n) exp 



.27rn ' 
~ % — ■ ( x ~ x ) 



(17) 



where n = (n x ,n y ,n z ) is 3-tuple of integers, /c 2 — (27r/L) 2 (n • n), and the term with n • n = is excluded from 
the summation. This is a direct consequence of substituting the known eigenmode functions of the Laplacian on T 3 , 
^„(x) = exp(z27rn • x/L), into eq. (||) for the correlation function. 

The method of images leads to the following alternative derivation of (|l7|). The correlation function on £ 3 is given 

by 

a(x,x') = -^3 J d 3 k P, & (fc)e ik ( x - x ') . (18) 

Assume, without loss of generality^ that x' lies in the Dirichlet domain with x as the basepoint. The contribution of 
each image has the form 

a(x,7ix') = J d 3 fcP$(fc) e lk '( x - x ') exp [i(k x ni + k y n{ + k z v? z )] . (19) 



4 For any two points, x and x' in A4 U , there exists a discrete motion 7 £ F such that the image 7X', would be in the Dirichlet 
domain with x as basepoint. Since x' and 7x' are equivalent for the compact space, considering the compact space correlation 
between x and 7x' instead of x and x' will give identical results. 
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from images 7,x' = x' + nPL. Summing over the contribution from the images, the correlation function on T 3 is 



1 - N N N 

v ; 7 n x =-N n v =-N n z =-N 



iknL 



(2tt) 3 



d 3 fcP$(fc)e ik(x - x,) lim 



sm 



[(AT + 1/2)^L] sin [(AT + l/2)fc„L] sin [(A/" + l/2)fc z L] 



sin(fc a; L/2) 



sin(fc !/ L/2) 



sin(fc z L/2) 



(20a) 



/OO 
d 3 fcP <I> (fc)e ik(x - x ' ) E <5(kL - 271-nj) 

= 73 Z^ p *(T n ) ex P 



,27rn ' 
-%— (x-x/) 



(20b) 



The final eq.(20b) is the same as (17) except that it contains a term P(n = 0) which is infinite for a wide class of 
power spectra; e.g.,P(k) cx fc Q , a < 0, including those a for which the integral J dkk 2 P(k) is convergent at k — 0, 
i.e., a > —3. 

The regularizing term can be easily calculated as well. Subtraction of the volume integral over N domains as in 
eq.fl) is described by the following substitution in eq. (20b) 



N N 



JY 



N N 



N 



iim y y y , 

n„= — N n„ = -7V n z = -N 



zknL 



lim [l-jo(fc*i/2)io(*„i/2)jo(feL/2)] E E E ' 

TV — ' oo * * ■ * 

n x = -N riv=-N n z =-N 



iknL 



(21) 



(where jo (a?) is the zeroth order spherical Bessel function) which leads to the following form for the regularized 
correlation function 



eT(x,x') = 7^)3/ d 3 fcP*(fcK k(x - x,) 



E 



<5(kL — 27rn,' 



24 



<y(k£) 



(22) 



As long as the power spectrum P$(fc) does not blow up faster that k~ 2 as k — > 0, the above regularization removes 
the zero mode contribution completely. In the case of the Harrison-Zeldovich spectrum (equal power per logarithm of 
k), where formally Pj>(fc) ~ fc~ 3 as k — > 0, the regularization suppresses the k — contribution but does not eliminate 
it completely. However, any physically motivated origin of P$(/c) such as from inflation does have an infrared cut-off, 
and all that is required of the regularization scheme is suppression of k = power. 

Ability to perform an analytic summation over all the images in T 3 would lead to the exact recovery of the positions 
of the discrete eigenvalues, 2n\n\/ L, and the eigenfunctions in this compact space. In more complex topologies (e.g., the 
CH spaces) one can only sum over a finite number of images and estimate the limit from that. If one tiles the universal 
cover out to A^ layers, one recovers the delta functions at k = 2im/L only approximately in the partial sum over 
images. The power of each discrete mode is aliased to a cubic cell of the reciprocal lattice^. In each cell, the kernel 
within square brackets in eq. (20a) has a peak at k — 2ir\n\/ L of height (2A^ + l) 3 and width w 1/AT, with damped 
oscillatory wings. Accurate regularization of the partial sum then consists of subtracting the total power in the 
reciprocal lattice cell around k = 0. 

The left panels of Figure ^ show the number density of states n T (k)/k 3 obtained by summation over N = 9 layers 
of images. The topmost panel shows the direct sum while the middle panel illustrates the effect of the regularizing 
term. The regularization procedure drastically reduces the spurious power below the fundamental frequency. The 
oscillatory wings of power aliasing inside each reciprocal cell can be eliminated if one averages over the results at each 
k obtained at different values of N. We found it is most effective to do it by Cesaro resummation: the effect is shown 
in the bottom panel. The result is a smoothed approximation to the underlying discrete spectrum (marked by arrows 
in the figure) with negligible contribution below the fundamental frequency and a positively defined spectrum. 

The plots of the right panel of Figure are analogous to those of the left panel except that they use the radial 
limiting procedure described in eq. (J24|) below, with = 9L. The radial limiting procedure does better at eliminating 



r.2w(m+^)/L 
'2n(m-h)/L 



dksin[(N + l/2)fcL]/ sin[fcL/2] = 2ty/L for any integer m and N. 
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the unwanted low k power at the resummation stage. The spectral lines are somewhat broader solely due to the fact 
that there are about tt/6 times fewer images within a sphere of radius r* = 9L than in the N = 9 layers. 

As we mentioned above, the regularization is not strictly complete for Harrison- Zeldovich like spectra, so the question 
of the numerically superior technique does arise. The resummation procedure effectively averages the sequence of 
partial sums up to a given distance. The superior low k power elimination of the radial limiting procedure is related 
to the fact that the product of the volume L 3 and the number of images within a radius r jitters around the volume 
47rr /3 on very short scales. The unwanted residual power in low k modes is distributed in this jitter, and is more 
readily removed by Cesaro resummation even if one does not go far in r. In the other case, the residual power in low 
k modes is distributed in a more orderly wave of wavelength ~ fc _1 in the sequence of partial sums. Hence, when one 
is summing images up to a finite distance d the residual power in modes with kd <C 1 is not averaged out. 




1 2 3 4 1 2 3 



wavenumber wavenumber 

FIG. 2. This illustrates the recovery of the known discrete spectrum in the T 3 model of size L = 2tt using the method of 
images. Plots of (2ir 2 /k 3 ) n T (k) in the three panels show the density of states obtained at different stages of our method. The 
left panels implement the limiting procedure of eq.(|^), varying the number of shells N of images that are summed. Images up 
to N = 9 layers (= 19 3 images ) are used. The right panels use the radial limiting procedure of eq.(^4j), with a cut-off radius 
r* that we use for CH spaces; here, r» = 9L, with only 3071 nearest images used. The vertical arrows mark the location of 
the discrete eigenvalues with the height proportional to the occupation number. The topmost panel shows the unregularized 
spectrum. Note the huge spurious contribution from the k — mode below the fundamental mode at k = 1. The middle panel 
shows the result of regularization. At this stage most of the spurious power is removed. The result of Cesaro resummation 
shown in the bottom panel demonstrates that our method finally recovers a smoothed power spectrum with negligible spurious 
long wavelength power. Note that the radial limiting procedure does better in terms of removing spurious power at low k. 

It is important to realize that at a given partial image sum the correlation function is obtained to far better accuracy 
than the power spectrum. The correlation function is a fc-space integral over the power spectrum and is accurately 
reproduced as long as power is sufficiently peaked around the correct eigenvalues and there is not much overlap 
between the adjacent peaks. The partial image sum generates a spectrum which can be approximated as the true 
spectrum convolved with a smearing function, W ((k — k t ) 2 / (Ak) 2 ) , around the true eigenvalues. Here (Ak) s quantifies 
the width of the smoothing. In a compact space no two points are physically separated by more than the diameter 
of the space. Hence, to get a reasonably accurate estimation of the correlation function, it is sufficient to ensure 
that R s ~ l/(Afc) s is larger than the diameter of the space. In the case of the simple T 3 model, in the Gaussian 
approximation to the smearing function W, the method of images with N layers gives R s — (2/ 3) dry/ N(N + 1) 
where dr — Lvo/2 is the diameter of the simple T 3 . Hence, the second layer (N — 2) approximation satisfies the 
R s > g?t condition by a comfortable margin, and even the first layer comes fairly close. 

The above estimate of the spread in the peaks around eigenvalues is based solely on the regularized spectrum done 
with summation in layers up to TV layers. We also carry out a Cesaro resummation to remove remaining spurious 
low k power. The resummation procedure effectively averages the cumulative results at each layer. It is not difficult 
to see that this reduces the amplitude of the peak by a factor of two and broadens the spectrum by the same factor. 
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Thus to reach the same spectral resolution after Cesaro resummation one need to go to N c « 2N layers. We prefer to 
implement the radial limiting procedure, which is simpler and does better at the resummation stage in suppressing 
residual low k power. Equating the number of images of the radial limiting procedure to that in N c layers implies 
choosing r*/L w (3/47r) s (2N c + 1). Hence, for the resummed spectrum obtained with the radial limiting procedure, 
the comfortable margin of N = 2 (which translates to iV c — 4) estimated above for obtaining accurate correlation 
functions translates to w 5.6L. In the second panel of Figure [|, we show the spectrum recovered at = 5.5L. 




12 3 4 

wavenumber 



FIG. 3. The recovery of the known discrete spectrum (2ty 2 /k 3 ) n T (k) in the T 3 model of size L = 2iv using the method 
of images. The top panel shows the regularized spectrum obtained at the lowest level approximation (N — 0, nearest image): 
only the low k cutoff is recovered. The second panel is the spectrum recovered after summing images up to a distance of 
r* = 5.5L, a level of approximation which should (and does) recover the correlation function quite well, although the spectrum 
is coarse. The bottom panel demonstrates the convergence of the regularized method of images: with a large enough number 
of images, the spectrum is also recovered with precision. In the lower two panels, Cesaro resummation has been carried out on 
the regularized spectrum. 

Another important remark is in order. In application to CMB, the positions of points x, x', between which the 
correlation function is to be calculated, are given by the length and directions of the photon path from the observer, 
i.e., correspond to the coordinates on the universal cover Ai u . In the method of images, the value of £J(x, x') for the 
pairs of points belonging to separate domains is found by symmetric replication of the correlation value computed 
after x' is mapped back into the Dirichlct domain around x. Thus, the method of images applied to the correlation 
function preserves at all levels of approximation the exact periodicity of the £J(x, x ') viewed as being defined on the 
universal cover. This periodicity is the major distinction between £|(x, x') an d the correlation function £$(x, x') in 
the simply connected universal cover space. In contrast, a correlation function determined as the inverse transform 
with respect to the universal cover eigenmodes of the approximate power spectrum will fail to obey the symmetries 
of the compact tiling strictly. This failure is greater the cruder the level of approximation used for computing the 
power spectrum. 

The success at reconstructing the correlation function by the method of images with even just a few layers of 



More precisely, when x' does not belong to the Dirichlet domain constructed around x. 
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summation is demonstrated in Figure |4|. In accordance with the estimate of the convergence of the method, summation 
over 5 layers (N=5.5) and above produces a correlation function almost indistinguishable from the exact result. 
TV = 3.5 is already very good and even N = 1.5 is an adequate approximation of the result. Even from the very 
beginning, with only the nearest image used, the method of images correctly catches the qualitative behavior of the 
correlation function in the compact space, which is dramatically different from the corresponding correlation function 
in non-compact flat space. 




45 90 135 180 

along the great circle on SLS 



FIG. 4. Correlation function in T 3 , calculated by the method of images with successively increasing number of images, is 
compared with the exact summation of eigenmodes of the Laplacian via Fourier transformation. Anticipating application to 
the CMB, we plot the correlation values between the points along the great circle of radius i? L s and one fixed point on the circle 
at (i? LS ,0, 0). The orientation of Cartesian coordinates coincides with directions of periodicity of the torus, so <\> is the polar 
angle. Monopole and dipole contributions along the circles (i.e., different on each circle) are subtracted. The results are shown 
for two values of R LS : The first case when R LS — 1 is a very symmetric one where points at 0, 90, 180, . . . degrees on the great 
circle are exact images of the point at the origin. The second case is a more general one. Labels mark the curves horn top to 
bottom at zero lag, <f> — 0. 



V. COMPACT HYPERBOLIC SPACES 

A. Correlation function on 7i 3 

The local isotropy and homogeneity of TC 3 implies £*(x, x') depends only on the proper distance, r = rf(x, x'), 
between the points x and x'. The eigenfunctions on the universal cover are of course well known for all homogeneous 
and isotropic models J3"l[ . Consequently £$(x, x') can be obtained through eq. (||). The role of £$(r) in the compact 
space calculation is to impose the desired power spectrum of $. In line with the application for which we developed 
our method, we now specialize to the scalar field <E> being the one describing cosmological gravitational potential 
fluctuations. 

The initial power spectrum of the gravitational potential P$(/c) is believed to be dictated by an early universe 
scenario for the generation of primordial perturbations. We assume that the initial perturbations are generated by 
quantum vacuum fluctuations during inflation. This leads to 
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where f3 = J (kd c ) 2 — 1, P<s>((3) = /3(/3 2 + l)P$(fc)/(27r 2 ) and, as before, r is in units of d c . In the simplest inflation 
models, the power per logarithmic interval of k, i.e.,V®, is approximately constant in the "subcurvature sector", 
defined by kd c > 1. This is the generalization of the Harrison- Zeldovich spectrum in spatially flat models to hyperbolic 
spaces |j32| , j33f . In appendix we outline a simplified calculation of the initial inflationary perturbation spectra V<& in 
hyperbolic geometry and derive a broader class of "tilted" spectra. Sub-horizon vacuum fluctuations during inflation 
are not expected to generate supercurvature modes, those with kd c < 1, which is why they are not included in eq. (p3|) . 
Indeed, since H 2 > l/(ad c ) 2 , for modes with kd c < 1 we always ha.vek/(aH) < 1 so inflation by itself does not provide 
a causal mechanism for their excitation. Moreover, the lowest non-zero eigenvalue, k\ > in compact spaces provides 
an infra-red cutoff in the spectrum which can be large enough in many CH spaces to exclude the supercurvature sector 



entirely (k\d c > 1). (See iVC.) Even if the space does support supercurvature modes, some physical mechanism 
needs to be invoked to excite them, e.g., as a byproduct of the creation of the compact space itself, but which could 
be accompanied by complex non-perturbative structure as well. To have quantitative predictions for P$(fc) would 
require addressing this possibility in a full quantum cosmological context. For a recent discussion of the creation of 
hyperbolic universes within a quantum cosmological framework see pi] . 



B. Numerical implementation of the method of images 

Equation (||) encodes the basic formula for calculating the correlation function using the method of images. For a 
numerical estimate, a limiting procedure such as eq. (Q) has to be used. In this form, the regularization term involves 
integrating £$(x, x') over a Dirichlet domain. Such terms can be numerically computed given the discrete group of 
motion, T, but it is usually cumbersome since the Dirichlet domains of CH spaces have complicated shapes which 
vary depending on the basepoint. 

A more effective and simpler radial limiting procedure in implementing the regularized method of images is to 
explicitly sum images up to a radius and regularize by subtracting the integral £|(r) over a spherical ball of finite 
radius r*: 



£J(x,x') = lim 



E 



4-7T 



dr sinh 2 r ££(r) 



= d(x,7j x ') ^ r j+i ■ 



(24) 



The volume element in the integral is the one for TC 3 . We have shown for the flat torus that this scheme works better 
numerically than formal regularization which subtracts integrals over Dirichlet domains. 

The plot on the left in Figure [5] illustrates the steps involved in implementing the regularized method of images. 
The value of £|, as a function of r* has some residual jitter, which arises because of the boundary effects due to the 
sharp "top- hat" averaging over a spherical ball chosen for the counterterm in eq. (p4|) . This can be smoothed out by 
resummation techniques |35j |. We use Cesaro resummation for this purpose. 

In hyperbolic spaces, the number of images within a radius r* grows exponentially for r*/ ' d c > 1 and it is not numer- 
ically feasible to extend direct summation to large values of r*/d c . The presence of the counterterm, however, besides 
regularizing, significantly improves convergence. This can be intuitively understood as follows: represents a 

sampling of a smooth function at discrete points rj. In a distant radial interval [r, r+dr], r i?>, dr ~ i?>, there are 
approximately (Att/Vm) sinh 2 r dr images. The sum, ^ r £$(rj), within this interval is similar to the (Monte-Carlo 
type) estimation of the integral, therefore one may approximate the sum over all distant images beyond a radius r* 
by an integral to obtain 



rj<r. 



47T 



dr sinh 2 r (r) 



(25) 



The tilde on £J(x, x') denotes the fact that it is approximate and unregularized. Subtracting the integral 
{Att/Ym) J °° dr sinh 2 r £$(r) as dictated by the regularization eq.(||), we recover the finite r* term of the limit- 
ing sequence in eq. (|24|) . This demonstrates that even at a finite r*, in addition to the explicit sum over images 



with rj < 



the expression for £|, in eq.(|24j) contains the gross contribution from all distant images with 



> r* 



Numerically we have found it suffices to evaluate the above expression up to r* 
i?> to obtain a convergent result for £$(x, x'). 



about 4 to 5 times the domain size 
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FIG. 5. The two plots illustrate the regularization of the correlation functions with and without an infrared cut-off at f3, > 0. 
The example shown is the correlation function at zero-point separation at some point on the CH manifold m004(-5,l). The 
plot on the left corresponds to £J,(x, x') where there is no f3 cut-off. The topmost panel shows the sampled values of £$(r) 
which contribute to the sum over images. The upper solid curve in the middle panel shows the divergent cumulative build up 
of the partial sum over images with successive addition of distant images. The dashed curve is the regularizing counterterm 
required to remove the zero-mode contribution and the lower solid curve is the cumulative value of the regularized partial sum 
£$(x, x'), which fluctuates around the true value once a sufficient number of images have been added. In the bottom panel 
this residual jitter (s°) in the estimate of the correlation is removed by Cesaro resummation. The sequence of lines, s 1 , s 2 and 
s shows the result of applying first, second and third order Cesaro resummation. The accuracy at the second order is usually 
sufficient. The plot on the right is analogous to the left but for the auxiliary correlation function £§(/?*, x, x') with an infrared 
cut-off at /3* = 4.0 below the first eigenvalue in the CH space. In contrast to the left panels, the cumulative sum over images 
is oscillatory which is more easily regularized, leading to much smaller residual jitter around the true value. 

Eq. (|2^) provides the simplest interpretation of our regularization procedure as an integral approximation to the 
total contribution of all distant images outside the region over which direct summation has been carried out. However, 
this interpretation may not be obvious in all cases. In fact, for the spectrum V{(3) = const that we use here, the 
correlation function on the universal cover Ti. 3 is given in terms of the hyperbolic sine and cosine integral functions, 
Shi(r) and Chi(r), respectively, by 

^(r) = Chi(r)-coth(r)Shi(r). (26) 

This is positive definite and does not fall off fast enough with r for its volume integral to converge. As a result, the 
integral in eq.(|2^) is not defined and the step from eq.(p4|) to eq.(|2^) is nontrivial, involving the regularization of an 
infinity which can be traced to the (3 = mode. To reinstate the intuitive interpretation in this case, we first show 
that it is valid for an auxiliary correlation function, r), which has an explicit infrared cutoff at /?» > 0. The 

auxiliary correlation function 

a(0 r)= r dp(3 MM vies 

**^*' ' Jp, (/? 2 + l) /3sinhr ' [P) 

= -Rc {[Ci [{i - /3„)r) + Ci ((* + /2 + * [Si [{i + /3„)r) + Si ((i - &)r)] /2 cothr} , (27) 

where Si(r) and Ci(r) are the sine and cosine integral functions, respectively. This function, £$(/?*, r), is no longer 
positive definite and its volume integral, although an improper integral, can be shown to be zero. 

The contribution to x, x ') °f the distant images with R ^> can now be evaluated explicitly: 

4tt r dr sinh\ a(P*,r) R ^ --^ sin (3*R + sin- 1 (k; 1 )) . (28) 

Jr Kdk„ R 

This is purely oscillatory with a zero mean. The plot on the right in Figure ^ shows that these oscillations of the 
regularization term precisely cancel out the oscillations (with growing amplitude) of the image sum, with no net effect 
on the limiting value of £|,(/3*,x, x'). 
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Having established that eq. ( [25| ) makes sense for £J(/3*, x, x'), where /3* > 0, the interpretation may now easily be 
extended to £|,(x, x ') by simply noting that given < ft* <C j3i, where j3\ is the wavenumber corresponding to the first 
eigenvalue of the Laplacian in the CH space, the method of images applied to the auxiliary r) an d £|(x, x') 

must converge to the same value. This result is demonstrated in Figure pi 
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p\=Q 




|S, = 5 



FIG. 6. Fhe analog of the bottom panels of Figure g for three values of the cut-off /3» . The more jagged curve is a regularized 
sequence of partial sums and the smoother curve is the result of a third order Cesaro resummation on the regularized result. 
The effectiveness of regularization is enhanced (reflected in the smaller amplitude of jitter) as one increases the value of j3* by 
directly cutting out more of the zero-mode contamination. The computation is robust and converges to the correct value as 
long as /3» <C Pi- For the CH manifold m004(-5,l) used in this figure, we estimate /3i > 5 (see Figure m. In the top three cases, 
/3» < 5, and the computed correlation converges to £ c (0) = 2 with high accuracy, independent of /3». When j3* > /3i, the value 
of £ c (0) is underestimated, as is seen in the bottom most panel where f3, = 8. 



C. Power spectrum 



In this section, we present our results on the power spectrum of CH manifolds obtained by applying the method 
of images to the density of states kernel n u (k,x, x') discussed in § III B . We emphasize that evaluating the density of 
states numerically in CH spaces is not our primary goal. What we present is a rough estimate of the power spectra 
in CH spaces that can be readily obtained as a 'byproduct' of our primary goal of computing correlation functions in 
CH spaces. 

Indeed, in the formalism of images, the singular delta-functions in n c (k) (see eq (11c)) are recovered, in principle, by 
the precise cancellation of the smooth contributions from all images. The volume of a sphere in 7i 3 , and consequently 
the number of images (of the Dirichlet domain of a CH space on the universal cover), grows exponentially with radius 
beyond r ~ d c . This is the primary constraint on the success of applying the method of images to the density of 
states. 

Our approximation for computing the CH space correlation function includes the gross integral estimate of the 
impact of distant images which results in a spread-out convolved density of states distribution. Increasing r* pro- 
gressively sharpens spectral profiles near the true positions of the discrete eigenvalues . Figure ^ shows a sample 
n c (k, x, x) for some of the CH manifolds at two random positions x in the first two columns. The third column shows 
the smoothed out estimation of the actual density of states, n c (k), for the CH manifolds that one can obtain by 
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averaging n c (/c,x, x) over ~ 10 3 points on the manifold. 

There is a definite signature of strong suppression of power at small (3 in all the cases that we have explored. 
This is qualitatively similar to the infrared cutoff known for the compact manifolds with flat and spherical topology. 
Quantitatively, the break appears around k ~ O (d^), consistent with the intuitive expectations. 




5 10 5 10 5 10 

p 

FIG. 7. The density of states exhibits suppression of long-wave power if the Universe has compact topology. Each row 
corresponds to one compact hyperbolic space, indicated by the label in the last column. Left and middle columns show 
n c (fc,x,x)/fc 3 computed with the method of images at two randomly selected points x. The last column shows the spatial 
averaged n c (k)/k 3 obtained by Monte-Carlo integration of n c (k, x, x)/fc 3 over the Dirichlet domain. The normalization is 
chosen such that the functions are equal to unity for the topologically trivial infinite open Universe. Vertical lines illustrate 
naive estimation of the cutoff at k ~ tt/cIm which holds surprisingly well for the cases checked. The first 6 examples of the CH 
space have about the same value of (1m- These spaces have volumes V/d 3 , around unity. The sixth space is m004(-5,l), the 
"small" space used in many of our examples. The seventh space is v3543(2,3), with volume V/n/d 3 = 6.45, used as an example 
of a "large" space in this paper. 



An infra-red cutoff at the lowest non-zero eigenvalue, k\ > 0, exists for all compact spaces. Cheeger's inequality |37 
provides a lower bound on k\ for a compact Riemannian manifold, M: 



k\ > — , 
- 2 ' 



he = inf 



A(S) 



s min{V{M 1 ),V(M 2 )} 



(29) 



where the infimum is taken over all possible surfaces, 5, that partition the space, M, into two subspaces, M\ and 
M2, i.e., M = Mi U M2 and S = dM\ = dM% (S is the boundary of M-y and M2). Cheeger's isoperimetric constant 
he depends more on the geometry than the topology of the space, with small values of he achieved for spaces having 
a "dumbbell- like" structure - a thin bottleneck which allows a partition of the space into two large volumes by a 
small-area surface p^ , ^8| . Regular shaped compact spaces do not allow eigenvalues which are too small. For example, 
the Cheeger limit for all flat T 3 manifolds is k\ > 2/L, where L is the longest side of the torus. Although the direct 
estimation of he is not simple, for any compact space M. with curvature bounded from below there exists a lower 
bound on he in terms of the diameter of M. p3j29[; for a 3-dimcnsional CH space, 



ki > h c /2 > — 



1/2 



eftcosli 0) 



0.92/4 



I'M- 



(30) 
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This result prohibits supercurvature modes for all CH spaces with cIm < 0.92 d c . 

There are other lower bounds on k\ that exist in the literature. In terms of the diameter alone, the bound 
k\ > 7r 2 /(2c? J v() 2 — max{— (n — l)K, 0} has been derived for any compact n dimensional space, with K = for flat 
geometry and K = ±d~ 2 for spherical and hyperbolic geometries, respectively f40j . For hyperbolic spaces, this bound 
is sharper than the one above for djw < 0.9 d c . There is another lower bound on k\ in terms of the volume as well 
as the diameter jy}; for CH spaces fci > VA</[27r<iM(sinh(2rfvi/d c ) — 2dvi/d c )]. This lower bound dominates in the 
supercurvature sector for volumes larger than w 7d^,. 




0.5 1 1.5 2 

diameter 



FIG. 8. Some known bounds on the first non-zero eigenvalue in CH spaces are shown. The heavily shaded region is the lower 
bound in eq. (^o|). The lightly shaded region is the lower bound derived in [40]. The dashed lines show the lower bound which 
depends on volume as well for three values of the volume. The shaded regions in the top right corner define two upper bounds 
on fci. The horizontally shaded region is the upper bound coming from comparison with Dirichlet eigenvalues [42]. The other 
region is the sharpest possible bound that can arise using the upper bound in terms of he [43] and the lower bound on he in 
eq. ©. 

There are also upper bounds on k\. The bound, k\ < 4hc/d c + 10lr^ pq ], does not allow for a firm conclusion 
that the space supports supercurvature modes, k\d c < 1, unless 4vi > 10.6 d c (using eq. (pp[)). Upper bounds based 
on comparison with the first Dirichlet eigenvalue ]42]| on a subdomain of M. cannot impose k\d c < 1, since Dirichlet 
eigenvalues cannot be less than d~ . 

The density of states n c (fc) defines the eigenvalue spectrum of the Laplacian on a compact space. It is well 
known that there exist very strong connections between the geometry and the topology of the space and the spectral 
properties of the Laplacian.^ As discussed above, the infra-red cut off in the spectrum is broadly determined by the 
size (diameter) of the space and its volume. 

The Weyl formula is an example of a general and powerful result connecting number of states to the topology of the 
compact space: for large k, the number of eigenvalues N(k) = #{j|fcj < k} up to a given value fc, in an n dimensional 



It is widely believed that, in two and three dimensions, the inverse problem of identifying a space from the spectrum of 
the Laplacian is spectrally rigid., i.e., there can be only a finite number of spaces which have the identical spectrum (modulo 
symmetries) [H. In higher dimensions, there are known exceptions, but spectral rigidity is quite generic. 
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compact space of volume Vm , is given by p9[ 

N(k) = #{j\kj <k}^ c(n)V M k n + O^™' 1 )) , c(n) = V(unit ball) / (2ir) n . (31) 

The constant c(n) is related to the volume within a unit sphere, set by the geometry of the space. In eq. (|3l|) , j counts 
the eigenvalues, including multiplicity. As a corollary, the eigenvalues asymptotically can be estimated: 

k, ~ ciny^U/VM) 1 ^ as j - oo . (32) 

The corrections to Weyl's asymptotic formula are £)(fc(™~ 2 ) +2 /(™ +1 )) for flat tori and 0(fc(™ _1 )/ ln(/c 2 )) for manifolds 
with negative curvature. 

Weyl's formula also points to the dependence of the infra-red cut off on size: the smaller the space, the larger the 
infra-red cutoff. Eq. (j32j) shows that in the limit of large eigenvalues the typical spacing between distinct eigenvalues 
is given by the inverse of its linear size. These are familiar and intuitive facts in Euclidean compact spaces such as the 
simple tori. Weyl's formula encourages the view that these broad features in the spectra of the compact manifolds 
transcend the geometry. Another important point of Weyl's formula relevant for our work is that the cumulative 
number of eigenstates asymptotically depends only on the volume of the manifold and not on its exact topology. 
Thus, the gross properties of the spectrum are shared by spaces with comparable volumes, and quantities which are 
fairly democratically weighted integrals of the power spectrum can be expected to be similar. 



VI. CORRELATIONS IN COMPACT HYPERBOLIC SPACES 

As a prelude to our computation of CMB temperature anisotropy correlations for CH spaces in fl6|| , we focus our 
attention in this section on the ^-correlation function ££(x, x') between points x and x' that belong to a 2-sphere 
around the origin in M. u — i.e., there exist discrete motions 7,7' S T such that the points 7X, 7'x' e A4 U are 
equidistant from the origin. 

Correlation functions of this kind arise in evaluating large angle CMB anisotropics associated with the gravitational 
potential on the sphere of last scattering (SLS) - the Naive Sachs- Wolfe (NSW) effect. (The radius of the SLS, 
i?LS ~ 2 tanh -1 yjl — Qq, is related to the density parameter, Slo-) The angular correlation C(q, q') between the NSW 
CMB anisotropy in two directions q and q' is given by 



C(q,q') = ^(q)^(q')^ = r LS )$fc, r M )> . 



(33) 



In simply connected universes, C(q,q') = C(q ■ q') is statistically isotropic. In contrast, for all compact universe 
models with Euclidean or Hyperbolic geometry, C(q, q') is statistically anisotropic. The breakdown of isotropy leads 
to characteristic patterns in the predicted CMB anisotropy determined by the shape of the Dirichlet domain around 
the observer. Moreover, except for the simple flat torus, the global inhomogeneity implies that the variance C(q,q) 
varies with direction q. This implies that the CMB sky would be a realization of an inhomogeneous field that would 
have characteristic patterns of 'loud' and 'quiet' regions. These two effects constitute two aspects of the signature of 
a compact universe. 

There are three distinct regimes for the correlation patterns on the sky in a compact space. For i?LS > R>, the 
patterns are mainly dominated by the mapping of the SLS into the compact space. The NSW-CMB patterns are 
characterized by spikes of positive correlation when the neighborhood of one point on the SLS is multiply-imaged on 
the SLS. For i?LS ^ R<, i.e., the SLS is well within a single domain, the compact space is indistinguishable from a 
simply-connected space with the same geometry. In the intermediate R < < i?LS ^5 R> regime, there is very little 
multiple imaging; nevertheless, the SLS is large enough to feci the compactness of the space. Typically the correlation 
pattern retains the structure of correlation in the simply connected space, but is significantly deformed. 
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FIG. 9. Typical angular correlation patterns C(qo, q) of the naive (or surface) Sachs-Wolfe effect on CMB anisotropy around 
a fixed direction go in the L(arge)CH model are shown as full-sky maps at the angular resolution of 5.2° x 5.2° pixels. The 
full-sky maps are plotted as pairs of 180° diameter hemispherical caps, one centered on the South Galactic Pole (SGP) and 
one on the North (NGP). go points to the NGP. In a simply connected universe, the contours of equal correlation would be 
concentric circles around go, due to the global isotropy. The three values of Qo are representative of the three regimes. In the 
top panel, the radius of the SLS, Rls, is greater than _R> and one sees multiple imaging. In the second, with i?< < Rls < R>, 
there are significant distortions but no multiple imaging. The bottom panel has R < -C Rls, and no observable correlation 
signatures of compactness. In all the maps, the dipole component of the correlation function has been subtracted. The maps 
have also been smoothed by a 1.66° Gaussian filter. 

Figure H illustrates the typical correlation pattern C(qo,q) in the large CH model v3543(5,l) around a fixed direction, 
<jo, which we chose to correspond to a fiducial "North Galactic Pole" (NGP). Of course, the global anisotropy implies 
that this pattern will differ for different choices of go- At Q,q = 0.6, the SLS is larger (R hs = 1.44) than the 
domain (i?> = 1.32), and spikes of enhanced correlation are seen with widely separated directions when the fixed 
point on the sphere, or points physically close to it, are multiply imaged on the SLS. At Qq = 0.8, the SLS is 
smaller (R < < R hs — 0.93 < i?>) and high correlation spikes due to multiple imaging are absent. Nevertheless, the 
compactness of the space is evident in the distorted contours around go- At Qq = 0.9, the SLS is completely contained 
within the domain, and the correlation around the NGP is circularly symmetric. As expected, in this regime the 
compactness of the space on scales much larger than the horizon has very little observational signal. The contours 
show only slight, observationally insignificant, distortions. 
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FIG. 10. The variations in the variance C(q,q) of a scalar field on three concentric spheres around the same basepoint. The 
radius of the spheres correspond to i?LS at the values of the density parameter Qo shown. The variance is shown in full-sky 
maps plotted as pairs of 180° diameter hemispherical caps, one centered on the South Galactic Pole (SGP) and one on the 
North (NGP). In addition to the COBE-DMRbeam smoothing, the maps have been smoothed by a 1.66° Gaussian filter. The 
loud feature seen at fio = 0.8 (R < < Rls < R>) corresponds to ~ 100% excess in the variance over its mean. At f2 = 0.6, the 
spot is multiply-imaged since the sphere is larger than the domain _Rls > R>- At Qo = 0.9 there is very little variation over 
the small sphere Rls <S R< ■ 

Figure |l0| plots the variance in the NSW-CMB temperature, C(q,q), in the large CH model v3543(5,l) for three 
values of Qq. The first two maps show a significant loud feature at f^o = 0.6 and fio = 0.8, corresponding to the radii 
i?LS = 1-44 and _Rls = 0.93, respectively. The loud spot on the sphere at flo = 0.8 is within the domain. At f2o — 0.6, 
the sphere is larger than the domain (i?> = 1.32), and the loud spot is multiply imaged on the sky. The third map has 
fio = 0.9, corresponding to a sphere of radius i?LS = 0.63, which is significantly smaller than the inradius, R < — 0.89. 
The variance does not show much variation over this small sphere around the observer. 

The CMB anisotropy has contributions other than the NSW term in eq. (|3^). In particular, there is an integrated 
Sachs- Wolfe contribution (ISW), with the integration along the line of sight from the SLS to the observer. The 
NSW term dominates when the density parameter f2o is close to unity. When the ISW contribution is important, it 
significantly modifies the correlation patterns, as discussed in our companion paper pa. 
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VII. DISCUSSION 



The computation of the angular correlation function of the CMB anisotropy requires machinery to compute spatial 
correlation functions on equal-time spatial 3-hypersurfaces. We have presented a general method of calculating these 
for multiply connected spatial sections which evades the task of eigenmode decomposition. This is particularly useful 
when considering compact hyperbolic models of the universe for which rather little is known about the spectrum of 
the Laplacian, and eigenmode decomposition is known to be difficult to obtain. 

We summarize here the basic knowledge we require of a given manifold to use our method and the steps to be 
followed to obtain a given accuracy level for the correlation functions. 

For the selected compact model Ai = M U JT of interest, we must be able to construct the tiling of the universal 
cover M. u from the generators of group T (§ ||) up to some distance r* from the origin, which will be determined by 
the required level of accuracy. This is in itself a daunting mathematical task, but, fortunately, for many thousand 
compact hyperbolic models, the Snappea package p8| gives enough information for us to carry out this step. 

Given the tiling, we perform our correlation function calculation in a sequence of radial shells of size r, te sting 
whether a stopping criterion based on a desired level of convergence in the regularized summation over images (§ VB) 
is satisfied; if so, this defines r*. It may be that the required is beyond the computational power at hand. For the 
correlation function calculations of interest for the CMB problem in the compact hyperbolic spaces we have tried, 10 5 
images are computationally very feasible (less than a day on a 433 MHz alpha workstation). This allows us to go out 
to pa 5d,M, more than adequate for convergence. 

It may be for manifolds with very many faces, the number of images required to converge could be prohibitive. 
Even if this is so, we obtain useful results out to the radius we can achieve because at each shell r the symmetries of 
the manifold are preserved in the correlation function. Indeed even the nearest few shells of images are enough for a 
qualitatively correct result on the basic pattern of correlations, which is also quantitatively not too bad. We showed 
explicitly for the flat torus (§ IV) that the correlation function is well determined as long as the density of states 
calculated for that number of shells has the power in each discrete mode localized to within Afc.; ~ 1/cIm around 
the true eigenvalue k t . Since the line of argument leading to this result is independent of the local geometry of the 
compact space, it applies to CH spaces as well. 

Our method has the merit of avoiding explicit eigenmode decomposition. Recently progress has been reported 
in directly computing the low-A: eigenmodes of the Laplacian for selected compact hyperbolic spaces 111]. It has 
proved necessary to use spectral line "deblend ing" t echniques in conjunction with these methods to get the eigenvalue 



spectrum accurately [|15| . As we showed in § |III B| , we can also use the method of images to calculate the spectrum 
of eigenvalues, though this is more difficult than correlation function evaluation since it requires longer summations 
to get narrow spectral lines, as is evident from Fig. [| and Q . These figures immediately suggesting deblending, but 
at this stage, only for the torus case do we explicitly know the shapes of the broadened spectral lines obtained by 
sequential image addition. One of the manifolds in which the low-A: eigenvalues were determined by Inoue fl4| is in 
common with those we have computed the spectrum for. We find his results agree with ours. 

The goals of this paper were to provide a detailed description of our regularized method of images, demonstrate 
that it works extremely well, and more generally show the basic effects of compactness on the correlation function 
of a scalar field and the density of states of the Laplacian. The computation of large angle CMB anisotropy, the 
features that arise in the CMB correlation that are characteristic of compact spaces, the detailed comparison with 
the COBE-DMRdata and the constraints that follow using large and small CH space examples are presented in the 
companion paper |l6) . A compilation of the constraints from the all sky COBE-dmr data on a large selection of CH 
spaces will be presented in |45fl. 
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APPENDIX A: INFLATIONARY PERTURBATIONS IN A HYPERBOLIC UNIVERSE 



It is widely accepted that the large scale structures in the present Universe grew out of small initial metric per- 
turbations via the mechanism of gravitational instability. The inflationary epoch in the Universe provides us with 
a setting for the generation of fluctuations with power on large scales. Cosmological perturbations are effectively 
massless (light) scalar fields residing on a background FRW space-time and the large scale power is related to the 
infra-red behavior of the light scalar field propagator in the inflationary epoch.^j In this Appendix, we present a 
simplified calculation of the initial power spectra for open inflation (hyperbolic geometry) to place the spectrum that 
we use in our work within a broad class of "tilted" inflationary spectra which are "scale free" on scales much smaller 
than d c and are analogous to the well-known tilted spectra in Euclidean (flat) inflationary universe. The scope of our 
analysis is limited by the particular choice of the vacuum state, as discussed below. 

Inflationary scenarios that lead to a simply connected hyperbolic universe generally require two stages JIG). This 
is because a single stage of inflation cannot provide homogeneity across the Hubble patch without inflating the local 
curvature in the patch to negligible values. A widely prescribed solution is to invoke a "creation stage" involving 
the creation of a universe with homogeneous hyperbolic spatial sections which is then followed by a sufficiently short 
inflationary stage (so as not to expand the curvature away) inflationary stage responsible for generating the primordial 
perturbations at cosmological scales. Since the second epoch of inflation cannot be long, the specifics of the first stage 
of creation of H 3 spatial sections could influence the spectrum of generation of perturbations at supercurvature scales. 
The possible effects on the spectrum have been extensively studied in recent literature jl7| in the context of open- 
bubble models in which the open universe resides in the interior of the bubble nucleated in the first-order phase 
transition from a larger inflating universe. It is doubtful, however, that the bubble mechanism can be used to produce 
compact hyperbolic spaces. More relevant is probably the mechanism of quantum creation of the universe "from 
nothing" where the topology of the created space is related to the topological properties of the instanton solution in 
the Euclidian section p4(| . Proper quantization of the fluctuations in such a scenario is, however, still a matter for 
the future. 

We assume that at the beginning of the second (inflationary) stage all quantum fields on the hyperbolic hypersurface 



are in the vacuum state and the vacuum modes in eq. (A6) are chosen by identifying the positive frequency part of 



the general solution to the mode evolution equation. This choice of the initial state sets the boundary of applicability 
for our considerations ^. It is important to bear in mind that, in general, t he creatio n st age could modify the proper 



choice of the vacuum state and consequently the predicted spectrum eqs. ( A7) and (jAjj). 

We outline the derivation of the power spectrum for a free scalar field ip which will be appropriately identified 
with metric perturbations subsequently. The scalar field equation allows the spatial and temporal dependence to 
be separated: ^(x, t) = V'fc( T ) G*k(x). The rapid expansion of the space-time during inflation would redshift away 
to insignificance the number density of any pre-inflationary quanta within the first few e-folds of expansion. As a 
result, it is a good approximation to assume that all the fields are in their vacuum state i.e., ipk(i~) = V'i + ' ) (' r )- The 
vacuum expectation value of the field, (iP 2 (t)}, is the coincidence limit (x' — » x) of the equal-time two-point function, 
(i/>(x', t)i/>(x, r)), and, hence, can be expressed as a mode-sum over the mode functions: 

(V> 2 } = ^3 J d 3 kMr) Mr) 

dk 
k 



^Mr)r k (r) . (AD 



The quantity within the square brackets represents the contribution per logarithmic fc— interval to the total power 
at a given time, r. We usually define the power spectrum by evaluating this quantity at some convenient time, r*, 
during inflation. For instance, t* could correspond to the time during the inflationary epoch when the comoving 
scale, fco, corresponding to the observed horizon was equal to the Hubble radius (&o|t:l| = fco T o)- The power spectrum, 
P^(k/ko), of the field, ij), is defined to be 



8 This infra-red behavior is linked to the infra-red divergence which strictly exists for massless scalars and for infinite duration 
of the accelerated phase. In the case of inflationary scenarios, both the finite duration and small effective mass serve to regulate 
the divergence. 

9 Although not explicitly used, the presence of the "creation stage" is essential. On the hyperbolic chart of the DeSitter space 
some observable 1-loop quantities computed in this vacuum would diverge on the (t — 0) boundary and hence the conditions 
assumed here at the onset of the inflationary epoch must breakdown in the past, presumably, due to the creation stage. 
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PAk/ko) = ^ Mn) Mn) ■ (A2) 

The present astrophysical comoving scales correspond to physical scales which were much larger than the Hubble 
radius at the end of inflation, hence it suffices to evaluate the amplitude of modes in the expression for the power 
spectrum, P^(fc/fco), in the limit, kr — *■ 0. 

We now compute the spectrum of initial perturbations generated in the second inflationary stage on a FRW model 
with H 3 spatial sections. 

The spatial modes of a scalar field on 7i 3 are given by orthonormal eigenfunctions of the Laplacian in spherical 
coordinates: |48| 



[3= V|k| 2 -1, r=|x|, .i = x/r, N t {0) = 



T{tf3 + £+l) 2 



T(if3+1) 



(A3) 



In the above equation, the P„(x) are the associated Legendre functions and T(v) is the Euler-Gamma function [(49[ . 
Owing to the isotropy of Tt 3 , the modes depend only k = |k|. Moreover, modes with kd c > f, (i.e., , (3 real and 
positive) form a complete orthonormal basis for free fields on 7i 3 , and it is convenient to use the wavenumber, (3, to 
label these modes. 

The temporal modes of the scalar field, tpp, on a hyperbolic universe, obey the equation 

^ + 2%' + [(p 2 + 1) + a 2 m 2 cS ] fa = , (A4) 

where m^ ff is the effective mass of the field and the prime denotes density derivatives with conformal time. We now 
consider the scalar field ip to be the fluctuations, S(f>, of the inflaton field |30 51 and then m^ ff is related to the 
quantities set by the evolution of the background inflaton field and the metric. 

The above setting allows us to investigate quantum fluctuations during inflation in some well defined regimes. The 
classification here closely follows that of inflation in a flat universe. We consider the following cases: 

i. uniform inflationary (deSitter) expansion and m 2 ff = 0. This is the analog of standard slow-roll inflation; 

ii. uniform inflationary expansion with a constant value of m 2 s /H 2 ^ 0. This is the analog of the class of models 
with inverted harmonic oscillator like potentials, such as natural inflation; and 

iii. constant equation of state during inflation. This is an abstraction of the power-law models of inflation. 

In each of the above, the conditions are meant to hold in an approximate fashion over the relevant range of astrophysical 
scales, just as for the flat inflation analogs. 

In the case of deSitter expansion, the scale factor is a = smh(Ht) / H = — cosechr/TJ, expressed as a function of 
cosmic time t and conformal time r. Eq. (A4) then reduces to 

^ - 2cothr^ + [(/3 2 + 1) + (m off /i7) 2 /sinhr] fa = 0. (A5) 
The general solution involves associated Legendre functions (coshr) and P ^ (coshr), multiplied by a 



factor sinh 2 r. The index /i = W9/4 — m 2 s /H 2 . The positive frequency (vacuum mode) solution is identified with 
the late time asymptotics ~ e~ tkT . The normalized vacuum mode solution, is 



^'-^ZVSi^-W^ <A6) 

where the overall normalization is determined by normalizing the probability current, i.e., setting the Wronskian of 
ip^ equal to ia~ 2 . 

T he s pectrum of initial perturbations in the inflaton field can be calculated from the vacuum mode solution through 
cq. (lAj): 
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Ah 



(A7) 



If Tn e ff — 0, the spectrum of initial perturbations in the infiaton field reduces to 



/3(/3 2 + l)- 



(A8) 



This is the analog of the Harrison-Zeldovich spectrum for hyperbolic models derived in ||32| , |33| . In the limit ,8 > 1, 
the spectrum goes over to the well-known flat "flicker noise" result. This is also the p ower spectrum which implies 
equal power per logarithmic interval in |k| that we use in our work on CH spaces (see § V A). 

The power spectrum in eq. (A.7) is best numerically computed by using Lancoz's formula for Gamma functions. 
Using the Sti rling approximation to Gamma functions for (3 3> 1, one recovers the standard tilted flat space spectrum 
from eq. dATh, "P* ~ fc 3 / 2 ^. 



0.9 
0.8 
0.7 

0.6 



Initial Power Spectra in Open — Inflation 
Tilted: /x= 1.7 

OPEN INFLATION 



-FLAT: % cxk 3 / 2 -^ 



OPEN: % oc£3/2- M 




FIG. 11. The initial power spectrum expected from inflation in a hyperbolic universe with a small tilt is shown as the solid 
curve (OPEN INFLATION). There is a sharp cutoff below kd c = 1 since the quantum fluctuations during inflation cannot 
excite supercurvature modes. The curve labeled FLAT shows the analogous tilted spectrum in a flat universe. The open 
inflation and the flat curves match at kd c ^> 1. Tilted initial spectra in Hyperbolic universes have often been taken to be a 
power law in kd c or /3. Although equivalent at kd c 3> 1, these differ at small kd c from the inflationary prediction; the former is 
labeled FLAT while the latter is depicted by the dashed curve labeled OPEN. The open inflation spectrum that we use in our 
paper corresponds to no-tilt, p = 3/2, is flat with a sharp cutoff at kd c — 1. 

In the more general case of constant equation of state, w = P/p, during inflation, the scale factor obeys a — 
sinh Q (r/a) where a = 2/(3w + 1). The deSitter expansion corresponds to the case w = — 1. The eq.(A4) for the 
temporal modes reduces to the form 



^ + 2 coth(r/a) ip'p + ((3 2 + 1)^ = 



(A9) 



The general solution of the above equation involves the same associated Legendre functions as the solution for equation 
( |A5| ) with index /i = 1/2 — a, rescaled wavenumber /3 2 — a 2 f3 2 = (A; 2 — l)a 2 , rescaled time f = r/oe, and multiplication 
by a factor of sinh M f . 
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The positive frequency (vacuum mode) solution is identified again with the late time asymptotics ~ e 
normalized vacuum mode solution, is 



-ikf 



The 



4 +) = h 2 



-=- smn r 



) ^ (coshf). 

-9+ 1 / 3 



(A10) 



The spectrum of pert urba tions is of the same form as (A7) with rescaled wavenumber j3 and ji = 1/2 — 2/(3w + 1). 
The spectra given by ( A7) are hyperbolic universe analogs of tilted spectra generated in flat universe inflation models 
such as power-law inflation. 



[1] J. R. Bond, D. Pogosyan and T. Souradeep, in: Proceedings of the XVIIIth Texas Symposium on Relativistic Astrophysics, 

ed. A. Olinto, J. Frieman, and D. N. Schramm, (World Scientific, 1997). 
[2] J. R. Bond, D. Pogosyan and T. Souradeep, Class. Quant. Grav. 15, 2671, (1998). 

[3] A. Dekel, D. Burstein and S. D. M. White, 1996, in Critical Dialogues in Cosmo logy, ed. N.Turok, World Scientific; D. N. 

Spergel, Class. Quant. Grav. 15, 2589, (1998); N. Bahcall and X. Fan, preprint ( |astro-ph/9804082| ). 
[4] G. F. R. Ellis, Gen. Rel. Grav. 2, 7,(1971). 

[5] D. D. Sokolov and V. F. Shvartsman, Zh. Eksp. Theor. Fiz. 66, 412, (1974) [JETP, 39, 196 (1974)] and references therein 

to the early papers on ghost searches. 
[6] M. Lachieze-Rey and J. -P. Luminet, Phys. Rep. 25, 136, (1995). 
[7] J. R. Gott, Mon. Not. R. Astr. Soc. 193, 153, (1980). 

[8] N.J. Cornish, D. N. Spergel and G. D. Starkman, Phys. Rev. Lett. 77, 215, (1996). 
[9] G. T. Horowitz and D. Marolf, J. High Energy Phys. 9807, 014, (1998). 
[10] G. Starkman, Class. Quantum Grav. 15, 2529, (1998); and other articles in the same issue: Proceedings of Topology and 

Cosmology, CWRU, Cleveland, Oct. 17-19, 1997 . 
[11] L. Z. Fang and H. Sato, Comm. Theor. Phys, (China) 2, 1055, (1983); H. V. Fagundes, Astrophys. J. 291, 450, (1985); 

ibid., Astrophys. J. 338 , 618, (1989); H. V. Fagundes and U. F. Wichoski, Astrophys. J. Lett. 322, L5, (1987). 
[12] H. V. Fagundes, Phys. Rev. Lett. 70, 1579, (1993); R. Lehoucq, M. Lachiese-Rey and J. -P. Luminet, Astron. Astrophys. 

313, 339, (1996); B. F. Roukema and A. C. Edge, Mon. Not. Roy. Astron. Soc. 292, 105, (1997); B. F. Roukema an d V. 

Blanloeil, Class. Quantum Grav. 15, 2645, (1998); J. -P. Luminet and B. F. Roukema, preprint, (astro-ph/9801364) and 

references therein. 

[13] R. Aurich and F. Steiner, P hysica D 39, 169, (1989); R. Aurich and J. Marklof, Physica D 92, 101, (1996); R. Aurich, 
preprint, ( |astro-ph/9903032| ). 

[14] K. Inoue, Class. Quant. Grav. 16, 3071, (1999). 

[15] N. J. Cornish and D. Spergel, preprint, ( |math.DG/9906017| ). 

[16] J. R. Bond, D. Pogosyan, and T. Souradeep in this volume (paper-II), (1999). 

[17] A. A. Starobinsky, JETP Lett. 57, 622, (1993); A. de Oliveira Costa, G. F. Smoot and A. A. Starobinsky, Astrophys. J. 
468, 457, (1996). 

[18] I. Y. Sokolov, JETP Lett. 57, 617, (1993); D. Stevens, D. Scott and J. Silk, Phys. Rev. Lett. 71, 20, (1993); A. de Oliveira 

Costa and G. F. Smoot, Astrophys. J. 448 477, (1995). 
[19] J. R.Bond, in Proc. of the XVIth Moriond Astrophysics Meeting, ed. F. R. Bouchet et al. ., (Editions Frontieres, France, 

1997); J. R. Bond, Proc. Natl. Acad. Sci. USA, 95, 35, (1998). 
[20] J. Levin, J. D. Barrow, E. F. Bunn and J. Silk, Phys. Rev. Lett. 79, 974, (1997). 
[21] D. D. Sokoloff and A. A. Starobinski, Sov. Astron. 19, 629, (1975). 
[22] J. A. Wolf, Space of Constant Curvature (5th ed.), (Publish or Perish, Inc., 1994). 
[23] E. B. Vinberg, Geometry II - Spaces of constant curvature, (Springer- Verlag, 1993). 

[24] W. P. Thurston, The Geometry of 3- Manifolds, lecture notes, Princeton University (1979); W.P. Thurston and J. R. Weeks, 

Sci. Am. 251, 108, (1984). 
[25] D. Gabai, R. G. Meyerhoff and N. Thurston, preprint MSRI 1996-058, (1996). 
[26] J. R. Weeks, Ph.D. Thesis, (Princeton, 1985). 

[27] S. V. Matveev and A. T. Fomenko, Usp. Mat. Nauk. 43, 5, (1988) [ Russ. Math. Surv. 43, 3, (1988)]. 

[28] J. R. Weeks, SnapPea: A c omputer program for creatin g and studying hyperbolic 3-manifolds, Univ. of Minnesota Geometry 

Center (freely available at http://www.geom.umn.edu). 
[29] I. Chavel, Eigenvalues in Riemannian geometry, (Academic Press, 1984). 



2G 



[30] G. Stark man, Class, Quant um. Grav. 15, 2529, (1998); J. Levin, E. Scannapieco, G. Gasperis, J. Silk and J. D. Barrow, 
preprint, [a.stro-ph/9807206; (1998); J. Levin, E. Scannapieco and J. Silk, Class. Quant. Grav. 15, 2689, (1998). 
E. Harrison, Rev. Mod. Phys. 39, 862, (1967). 
D. H. Lyth and E. D. Stewart, Phys. Lett. B252, 336, (1990). 
B. Ratra and P. J. E. Peebles, Astrophys. J. Lett. 432, L5, (1994). 

G. W. Gibbons, Nucl. Phys. b f B 472, 683, (1 996); G. W. Gibbons, Class. Quantum Grav. 15, 2605, (1998); S. S. e Costa 
and H. V. Fagundes, preprint, |gr-qc/9801066; , (1998). 
G. H. Hardy, Divergent Series, (Oxford University Press, London, 1956). 
N. L. Balazs and A. Voros, Phys. Rep. 143, 109, (1986). 

J. Cheeger, in Problems in Analysis, (A Symposium in honor of S. Bochner), (Princeton University Press, 1970). 
P. Buser, Proc. Symp. Pure Math. 36, 29, (1980). 

P. H. Berard, Spectral Geometry: Direct and Inverse Problems, Lec. Notes in Mathematics, 1207, (Springer- Verlag, Berlin, 
1980). 

P. Li and S. T. Yau, Proc. Symp. Pure Math. 36, 205, (1980). 
S. T. Yau, Ann. scient. Ec. Norm. Sup. 8, 487, (1975). 
S. Y. Cheng, Math. Z. 143, 289, (1975). 
P. Buser, Ann. scient. Ec. Norm. Sup. 15, 213, (1982). 

D. Gurarie, in Forty More Years of Ramifications: Spectral Asymptotics, ed. by S. A. Fulling and F.J. Narcowich, (Texas 
A&M University, College Station, Texas, 1991). 
J. R. Bond, D. Pogosyan and T. Souradeep, in preparation. 

J. R. Gott, Nature, 295, 30, (1982); A. H. Guth, and E. J. Weinberg, Nucl. Phys. B, 212, 321, (1983); M. Bucher and N. 
Turok, Phys. Rev., D 52, 5538, (1995). 

M. Sasaki, T. Tanaka, K. Yamamoto, and J. Yokoyama, Prog. Theor. Phys. 90, 1019, (1993); K. Yamamoto, M. Sasaki, 
and T. Tanaka, Astrophys. J. 455, 412, (1995); M. Bucher, A. S. Goldhaber and N. Turok, Phys. Rev. D52, 3314, (1995); 
K. Gorski, et al. , Astrophys J. Suppl., 114, 1, (1998). 
L.F. Abbott, and R. K. Schaeffer, Astrophys. J. 308, 546, (1986). 

M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, (Dover, New York, 1970). 
V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rep. 215, 203, (1992). 

J. R. Bond, Theory and Observations of the Cosmic Background Radiation, in Cosmology and Large Scale Structure, Les 
Houches Session LX, August 1993, ed. R. Schaeffer, (Elsevier Science Press, 1996). 



27 



